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AUTOMATIC COLOR BALANCE 
BACKGROUND OF HE INVENTION 

1 , Field of the Invention 

5 The present invention relates to the field of imaging, the field of computer assisted 

imaging, the field of digital imaging, and the field of automatically controlled 
enhancement of specific attributes of digital imaging data such as contrast. 

2. Background of the Invention 

10 The proliferation of digital imaging means such as photography, scanning, copying, 

printing and digital cameras has resulted in a large volume of color imagery. Since none 
of the devices produce consistently perfect color, especially in the hands of unskilled 
amateurs, there is a need to correct the color of images. Color correction has been the 
object of much effort in the imaging art but there remains a need for simple correction 

15 methods that can be applied in an automated way. 

One approach to correcting color is based on the idea that improper colors in the scene 
are the result of illumination that is not perfectly white. If the illuminant of the scene can 
be determined, then the colors of the scene can be corrected to their appearance under 

20 some standard reference illuminant. Many methods make the assumption that specular 
reflections in the image have a color corresponding to the illuminant so that the brightest 
pixels in the image can be used to recover the illuminant color. Examples of this 
approach include the following. U.S. Patent No. 4,685,071 describes estimating the 
illuminant by determining the locus of intersection of lines fitted through sets of points of 

25 constant hue and varying saturation. US Pat. 5,495,428 teaches a similar method in which 
the improvement involves weighting the lines according to their reliability. U.S. Patent 
No. 5,825,916 involves a related fitting of lines to a smoothed chromaticity bitmap. U.S. 
Patent No. 6,104,830 discloses a similar procedure in which the location of the lines is 
estimated by means of a Hough transform. In "Signal Processing by the Input Interface to 

30 a Digital Color Laser Copier" by A. Usami, SID 90 Digest, p. 498-500 ( 1 990) the 

brightest non-white pixel is considered representative of the illuminant and is adjusted to 
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a neutral color. However, these methods fail when detector saturation occurs since the 
brightest pixels then no longer represent the illuminant and, except for the Usami 
procedure, are computationally expensive. 

5 An alternative approach to illuminant estimation is to examine the gamut of colors in an 
image. If certain colors are present in the image particular illuminants can be excluded. 
For example, since objects are usually colored by virtue of reflecting light, if a scene 
contains the color red then the illuminant must contain red and cannot, for instance, be m 
blue. The gamut mapping procedure for illuminant recovery is described in "A novel 

10 algorithm for color constancy' 5 by D. Forsyth, Int. J. Comput. Vision, 5, p. 5-36 (1990). A 
more efficient version in a chromaticity color space has been developed by G. Finlayson 
and S. Hordley, Proc. IEEE Conf Comput. Vision Pali. Recogn., p. 60-65 (1998). 
European Pat. 0 862,336 teaches the use of the method in a digital camera. These 
methods are computationally intensive an do not uniquely identify the illuminant with 

15 additional assumptions. Moreover, digital images can be subject to arbitrary color 

manipulation so that color imbalance does not necessarily result from illuminant changes. 

Yet another method of color correction is based on the gray world assumption introduced 
by Evans in U.S. Patent No. 2,571,697. The method relies on the idea that in a complex 

20 natural scene such as typically occurs amateur photographs the average of all the colors is 
gray. Thus, by adjusting the mean color of the image to gray, color correction can be 
achieved. However, this method fails when the scene content does not, in fact, 
correspond to an average gray. This happens, for instance, when an object of a single 
color dominates the scene or in digital images, such as business graphics, which have a - 

25 simple color distribution. There have been attempts to improve the method by applying it 
to individual luminance ranges in the image. Examples include U.S. Patent No. 
5,233,413, U.S. Patent No. 5,357,352 and U.S. Patent No. 5,420,704. Another variation 
in U.S. Patent No. 5,926,291 seeks to use only colors of low and high brightness and also 
low saturation as the basis for correction. The same gray world assumption is used in the 

30 retinex family of algorithms as discussed in "Investigations into multiscale retinex", K. 
Barnard and B. Funt, Color Imaging in Multimedia '98, p. 9-17, Derby, UK, March 1998. 
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None of these methods are, however, fully satisfactory because of failure of the gray 
world assumption. 

In order to improve color balancing performance some workers have taken advantage of 
5 the statistical distribution of image types submitted to the color correction system and 
have developed corrections tailored to certain common types of color defects. Examples 
include U.S. Patent No. 4,339,5 1 7 and U.S. Patent No. 6^097,836. However, such 
approaches are useless when the images to be processed do not fall into a few simple 
categories. Other correction methods attempt to capture the experience of imaging 
10 experts by developing rules for image correction based on examination of a very large 
number of images as exemplified by U.S. Patent No. 5,694,484. In WO 97/01 151 there is 
disclosed a color correction system that is taught color preferences through the use of 
neural networks. Such methods frequently fail when they encounter images not in the 
original training set. Moreover, the effort of developing such a method is very great 
15 because of the large number of images on which it is based and the resulting method is 
hard to understand and modify because of its complexity. In the case of neural networks, 
there is a danger of over-training, where correction of the training set improves at the. 
expense of generality in the correction performance. 

20 A range of empirical color correction methods also exist, which are based on statistical 
analysis of color histograms and sometimes also brightness histograms. Examples include 
U.S. Patent No. 4,729,016, U.S. Patent No. 4,984,071, U.S. Patent No. 5,1 17,293, U.S. 
Patent No. 5,323,241 and U.S. Patent No. 6,151,410. Most of these methods place 
particular emphasis on the highlight and shadow regions of the histogram though some, 

25 such as U.S. Patent No. 6,151,410, specifically exclude some of these regions on the 

grounds that the data are unreliable. These methods depend on the image data set used to 
derive the statistical analysis of the histogram and can be unsatisfactory for images of a 
type not in the original data set. At least some cause over-correction of the image when 
applied repeatedly and can result in information loss through clipping of the lowest and 

30 highest intensities in the image. These methods are, further, inherently incapable of 
correcting for different scene illuminants. 
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There are also color correction methods based on manually specifying a black or white 
point in the image or a color that should be considered as neutral gray. Such a capability 
is available as software in the "Balance to sample" feature of PhotoStyler 2.0 (Aldus 
Corporation, 41 1 First Avenue South, Seattle, WA 98104), in the "Curves" feature of 

5 Photoshop 5.5 (Adobe Systems Incorporated, 345 Park Avenue, San Jose, California 
951 10-2704) and in the "Automatic" mode of the "Tint" feature in PhotoDraw 2000 
(Microsoft Corporation, One Microsoft Way, Redmond, WA 98052-6399). Correction 
using manually specified highlight and shadow regions is disclosed in U.S. Patent No. 
5,062,058 and the falls within the claims of U.S. Patent No. 5,487,020. These methods, 

10 however, require manual intervention and it can be difficult to select the optimal black 
and white points to achieve the desired correction. There have been attempts to automate 
this correction process. In 1996 Photoshop 4.0 introduced the "Auto Levels" feature that 
stretches the three color channel histograms to full range, by default clipping the top and 
bottom 0.5% of the channel values. A similar feature is available as "Auto Tonal 

15 Adjustment" in PhotoStyler 2.0. Additionally, U.S. Patent No. 5,812,286 teaches such a 
method of correction. These methods have the disadvantage that part of the image 
information is lost through the clipping process and, further, image contrast is 
undesirably changed along with the correction of color. An attempt to solve this difficulty 
is disclosed in U.S. Patent No. 5,371,615 wherein the RGB color triplets for each image 

20 pixel are examined to determine the blackest non-black pixel as min[max(R,G,B)] and 
the whitest non-white pixel as max[min(R,G,B)] ignoring exactly black and exactly white 
pixels. Subsequently a black point'W m j n is specified as having all three color channels 
equal to min[max(R,G,B)] and a white point W max as having all three color channels 
equal to max[min(R,G,B)] an d then each color channel value Xj n is corrected to a new 

25 value X out according to: 

X 0 ut = (W max — Wmin) x Q^'m — X m j n ) / (X max — X m j n ) + W m j n . 

This procedure, however, has the disadvantage that the color correction can depend on as 
few as two pixels in the image. This renders the method susceptible to noise and to 
defective pixels in digital camera detectors. At the same time the method does not retain 
30 the contrast of the image, an effect that can be especially marked when the green channel 
does not participate in the definition of W min and W max . 
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SUMMARY OF THE INVENTION 

It is an aspect of this invention to provide a method of color correction for digital images. 
A further aspect of the invention is to provide a method of removing strong color casts 
from images. Another aspect of the invention is to provide a means of automatically 
5 correcting color balance in images. Yet another aspect of the invention is to correct color 
balance in an image by a means that allows compensation for the color of the scene 
illuminant. Still another aspect of the invention is to provide a method of color correction 
that leaves image contrast substantially unchanged. 

10 These and other aspects of the invention are achieved by averaging color channels in 
regions near the minimum of the histograms of the color channels and selecting the 
smallest of the average color values as a black point; averaging color channels in regions 
near the maximum of the histograms of the color channels and selecting the largest of the 
average color values as a white point; correcting the color channels by adjusting the 

15 aforementioned color averages to the match the values of the black point and white point, 
this adjustment being done optionally in a way that limits clipping of color values; . 
optionally providing an illuminant correction; and optionally providing a method to 
restore image contrast after color correction. 



Figure 1 shows brightness histograms with identified regions of maximum color values 
and minimum color values, average points, black points, white points, and color value 
limits. 

25 Figure 2 shows brightness histograms with identified regions of maximum color values 
and minimum color values, average points, black points, white points, and color value 



30 The invention operates on digital images in which information is supplied as a grid of 
points or pixels. Each pixel is associated with several numbers representing color 
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information. Typically, these are three 8- or 12-bit values representing respectively red 
(R), green (G) and blue (B). While the invention will be illustrated using such RGB 
images as examples, it will be understood that more than three or more colors can be used 
to represent color information, that the color space used to define color need not be 
5 restricted to RGB and that the color information is not restricted to any specific form of 
numeric representation. Although the invention will be described in terms of calculations 
using all the image pixels, it also understood that the invention may be applied to only a 
portion of the image pixels, either with respect to the analysis used to derive the means of 
correction or with respect to the regions of the image that undergo the correction, or both. 
10 Whereas the invention will be described in terms of certain illustrative embodiments it 
will be understood that wide variation of specific realizations of the invention is possible 
within its scope. 

In the practice of the invention, there are four sub-component processes that may 
15 combined, each of the sub-components of the process being novel. The first two sub- 
component processes relate to an initial color correction step in data format, and the 
second two steps effect actual implementations of the color adjustments on the data or 
images. The initial color correction, in general terms involves averaging at least two 
color channels in regions near the minimum of histograms of the at least two color 
20 channels; selecting the of the average color values as a black point or dark reference 
point; averaging at least two color channels in regions near the maximum of the 
histograms of the at least two color channels and selecting the largest of the average color 
values as a white point or light reference point; correcting the at least two color channels 
by adjusting the smallest and the largest color averages to the respectively match the 
25 values of the black point and white point to form corrected image data. It is preferred to 
impose a clipping limit restriction on the process of matching average color values to the 
reference points (i.e., the black point and white point) and the first two sub-component 
processes for initial color correction differ in how this may be accomplished. 
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The first inventive step will be referred to as the "direct color correction." The process 
would ordinarily be performed for each color (e.g., R, G and B), but may be performed 
for only two colors and provide an effect on an image. The process may comprise: 

1) For each color, forming a histogram of the darkest pixels in the image data 

5 (ignoring pixels with 0, 0, 0 values, that is true black pixels). The measure of 

darkness is the largest of the color values at a pixel. 

2) For each color, forming a histogram of the lightest pixels in the image data 
(ignoring pixels with 255, 255, 255 values, that is true white pixels). The 
measure of lightness is the smallest of the color values at a pixel. 

10 3) For the histogram of darkest pixels for each color determining a number (that 

is, more than one, or at least two) of the darkest pixels and compute a dark 
average value for this number of pixels. It is preferred that the number be a 
specific fraction of the area under the histogram. 

4) For the histogram of lightest pixels for each color determining a number (that 
15 is, more than one, or at least two) of the lightest pixels and compute a light 

average value for this number of pixels. It is preferred that the number be a 
specific fraction of the area under the histogram. 

5) Selecting the smallest dark average from among those determined for each 
channel for the darkest pixels as a dark reference point. 

20 6) Selecting the largest light average from among those determined for each 

channel for the lightest pixels as a light reference point. 

7) Specifying a fraction of image pixels as a maximum clipping fraction for each 
color channel for dark colors and for light colors. This defines a dark clip limit 
value for each color and a light clip limit value for each color. Any color value 

25 less than the dark clip limit value may be set to the minimum possible value 

for the color (e.g., 0). Any color value greater than the light clip limit value 
may be set to the maximum possible value for the color (e.g., 255). 

8) Conceptually, for each channel displacing the dark average value towards the 
dark reference (where one of these averages is the dark reference and does not 

30 move). The movement continues until either the average aligns with the 
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reference or the dark clip limit value aligns with the minimum channel value, 
whichever comes first. 

9) Conceptually, for each channel displacing the light average value towards the 
light reference (where one of these averages is the light reference and does not 
move). The movement continues until either the average aligns with the 
reference or the light clip limit value aligns with the maximum channel value, 
whichever comes first. 

10) For all the values of each color that are not clipped, constructing a linear 
look-up table between a point defined by the original and displaced positions 
of the dark average and a point defined by the original and displaced positions 
of the light average and applying the look-up table to modify the color values 
in the image. 

Elements of this process are diagramed in Figure 1. The top panel shows the original 
histograms for red, green and blue channels. The shaded areas at each end of the 
histogram represent the fraction of pixels used to compute the dark and light averages. 
The values of these averages are marked with short solid vertical lines. The dark 
reference is marked DR and the light reference is marked LR. The center panel shows 
perfect alignment of all three dark averages with the dark reference and of all three light 
averages with the light reference, as indicated by the dotted lines. Under these conditions 
low values of the red channel would be excessively clipped (EC) to a degree greater than 
permitted by the specified maximum clipping fraction. Similarly high values of green 
channel would be excessively clipped (EC). The bottom panel shows the results of the 
above specified process, where the dark average or the red channel is moved only as far 
towards the dark reference as is consistent with the permitted amount of clipping PC. 
Thus, the low averages of the green and blue channels align with the dark reference but 
that of the red channel does not. For the red channel the dark clip limit value aligns with 
the minimum color value. Similarly, the light average of the green channel is moved only 
as far towards the light reference as is consistent with the permitted amount of clipping 
PC. Thus while the light averages of the red and blue channels align with the light 
reference, for the green channel the light clip limit aligns with the maximum color value. 



The second inventive step will be referred to as the "color cast correction." The process 
would ordinarily be performed for each color (e.g., R, G and B), but may be performed 
for only two colors and provide an effect on an image. The process may comprise: 

1) For each color, forming a histogram of the darkest pixels in the image data 
(ignoring pixels with 0, 0, 0 values, that is true black pixels). The measure of 
darkness is the largest of the color values at a pixel. 

2) For each color, forming a histogram of the lightest pixels in the image data 
(ignoring pixels with 255, 255, 255 values, that is true white pixels). The 
measure of lightness is the smallest of the color values at a pixel. 

3) For the histogram of darkest pixels for each color determining a number (that 
is, more than one, or at least two) of the darkest pixels and compute a dark 
average value for this number of pixels. It is preferred that the number be a 
specific fraction of the area under the histogram. 

4) For the histogram of lightest pixels for each color determining a number (that 
is, more than one, or at least two) of the lightest pixels and compute a light 
average value for this number of pixels. It is preferred that the number be a 
specific fraction of the area under the histogram, 

5) Selecting the smallest dark average from among those determined for each 
channel for the darkest pixels as a dark reference point. 

6) Selecting the largest light average from among those determined for each 
channel for the lightest pixels as a light reference point. 

7) Specifying a fraction of image pixels as a maximum clipping fraction for each 
color channel for dark colors and for light colors. This defines a dark clip limit 
value for each color. and a light clip limit value for each color. Any color value- 
less than the dark clip limit value may be set to the minimum possible value 
for the color (e.g., 0). Any color value greater than the light clip limit value 
may be set to the maximum possible value for the color (e.g., 255). 



8) Conceptually, for each channel displacing the dark average value towards the 
dark reference (where one of these averages is the dark reference and does not 
move). The movement continues until there is perfect alignment. 

9) Conceptually, for each channel displacing the light average value towards the 
light reference (where one of these averages is the light reference and does not 
move). The movement continues until there is perfect alignment. 

10) Conceptually, and if required, displacing the aligned dark averages in unison 
to higher values until the low values of every color are clipped to an extent no 
greater than the maximum clipping fraction. 

1 1) Conceptually, and if required, displacing the aligned light averages in unison 
to lower values until the high values of every color are clipped to an extent no 
greater than the maximum clipping fraction. 

12) For all the values of each color that are not clipped, constructing a linear 
look-up table between a point defined by the original and displaced positions 
of the dark average and a point defined by the original and displaced positions 
of the light average and applying the look-up table to modify the color values 
in the image. 

Elements of this process are diagramed in Figure 2. The top panel shows the original 
histograms for red, green and blue channels. The shaded areas at each end of the 
histogram represent the fraction of pixels used to compute the dark and light averages. 
The values of these averages are marked with short solid vertical lines. The dark 
reference is marked DR and the light reference is marked LR. The center panel shows 
perfect alignment of all three dark averages with the dark reference and of all three light 
averages with the light reference, as indicated by the dotted lines. Under these conditions 
low values of the red channel would be excessively clipped (EC) to a degree greater than 
permitted by the specified maximum clipping fraction. Similarly high values of green 
channel would be excessively clipped (EC). The bottom panel shows the results of the 
above specified process, where the dark averages maintain alignment but are displaced to 
higher values to reduce the clipping of the red channel to the permitted amount of 
clipping (PC). Similarly the light averages maintain alignment but are displaced to lower 
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values until the green channel has no more than the permitted amount of clipping (PC). 
Thus while all the averages are aligned in the center and bottom panels, the separation 
between the dark and light averages or } equivalently the separation between dark and 
light references, is reduced in the bottom panel. 

5 

This process can be alternatively and more mathematically described as the pixels of the 
image being examined with regard to the RGB values at each pixel in order to define 
highlight and shadow regions of color. This may be done in a number of ways, for 
example by building histograms of color distribution in each channel. If these histograms 

10 are computed in cumulative (i.e., integral) form it is possible to define a shadow region as 
those colors lying below some integral threshold and highlights as those colors lying 
above some integral threshold. Such thresholds correspond to fractions of the total 
number of pixels in the image. It is also possible to use only a single channel, for instance 
a brightness channel, to define highlight and shadow regions. In practice a convenient 

15 way to define the shadow and highlight regions is the following. A histogram is 

constructed of the minimum values, min[R(i,j), G(i,j), B(i j)], in each color triplet for 
every pixel with coordinates (i j) and pixel color values belonging to a certain integral 
fraction measured from the high end of the histogram are defined as belonging to the 
highlight region. Here min[R(i,j) } G(i,j), B(i,j)] means the smallest value of the three 

20 color values defining a pixel and the aforementioned fraction is termed a region fraction. 
Similarly, a histogram is constructed of the maximum values, max[R(i,j), G(i j), B(i j)], in 
each color triplet for every pixel and pixel color values belonging to a certain integral 
fraction measured from the low end of the histogram are defined as belonging to the 
shadow region. Here max[R(i,j), G(i,j), B(i,j)] means the largest value of the three color 

25 values defining a pixel. Though pure black pixels and pure white pixels can be included 
in these calculations, it is preferred to omit them since these pixels are often unreliable 
because of the prior processing history of the image. In the event that the histograms are 
empty because the image contains only black, or only white, or only black and white 
pixels, correction is not performed. The region fractions used to define shadows and 

30 highlights can be the same or different, though it is convenient to use the same fraction 
for both ranges. The region fraction can be a number greater than 0% and less than 50%, 
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with a preferred range of 0.04% to 40% and a most preferred range of about 4%. It has 
been found that a range of about 4% is suitable for balancing the color of the majority of 
images encountered, without any user intervention, in an automatic way. Very high 
values near the top of the range can be used for extremely poor images with exceptionally 
5 poor color balance. For this reason, though the invention can function without user 
intervention for normally encountered images, it is useful to also provide the user with 
control of the fraction defining shadow and highlight regions. The region fraction could 
be set directly by the user but it has been found that a non-linear control best represents 
the visual effect of changing this fraction. For example, an implementation where linear 
10 steps in the control correspond to logarithmic steps in the region fraction has been found 
effective. 

A description of the second step would be that, for all pixels whose max[R(i,j), G(i j), 
B(i j)] values lie in the shadow region, the individual color channels are averaged to 

15 produce three average colors Ri ow , G| 0W and Bi ow constituting histogram control points. 
The minimum of these three values, Xi ow , is then computed as a reference point. 
Similarly, for all pixels whose min[R(i j), G(i,j), B(i,j)] values lie in the highlight region, 
the individual color channels are averaged to produce three average colors R h i g h, G h i g h, 
Bhigh constituting additional histogram control points. The maximum of these three 

20 values, Xhigh, is then computed as a second reference point. The intent of the invention is 
to stretch the color histograms in such a way that the control points R| 0W , G] 0W and Bi ow 
coincide with the Xi ow reference point and the control points Rhigh, G h i g h, B h i g h coincide 
with Xhigh reference point, though this intent may be restricted by conditions to be 
described later. This histogram stretching can be conveniently and efficiently achieved by 

25 construction of look-up tables that relate initial channel value to channel values after 
correction. The use of color averaging confers noise immunity on the algorithm. 

As will be discussed later in the description of specific embodiments of the invention, the 
alignment of histogram control points with their respective reference points is restricted 
30 by control of histogram clipping. By clipping is meant the conversion of a range of 
shadow or highlight input values into a single shadow or highlight output value. This 
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results in information loss from the image and is consequently undesirable. A clipping 
fraction can be defined as that range of color values which are responsible for a particular 
fraction of the of the area under the histogram, which in turn corresponds to a particular 
fraction of pixels in the image. The clipping fraction is measured from the ends of the 
5 color channel histograms in a manner similar to that described for the region fraction. It is 
a measure of the amount of information that is lost from the image as a result of clipping. 
The clipping fraction can be different at the high and low ends of the histogram but it is 
preferred to use the same fraction at both ends. The clipping fraction can be a number 
from 0% to less than the region fraction, with a preferred range of from about 0.0001% to 
10 about 10%, and a most preferred range of from about 0.01% to about 5%. Additionally, it 
is especially preferred to make the clipping fraction dependent on the region fraction such 
that the clipping fraction is from about 0.02 to about 0.2 times the region fraction. Most 
especially preferred is a clipping fraction of about 0.2 times the region fraction. For the 
preferred region fraction of about 5% the preferred clipping fraction is thus about 0.1%. 

15 

First Embodiment 

This embodiment is suitable for processing the majority of normal images that required 
color balancing, especially when it is desired to perform the correction automatically 
without user intervention. In the shadow region each of the control points Ri 0 w, Gi 0W) Bi ow 

20 that lies above Xi ow is moved towards the reference point Xi ow , subject to the condition 
that the clipping fraction is not exceeded. Thus, if a particular channel does not exceed 
the predetermined clipping fraction its control point is aligned with Xi ow . Otherwise, the 
channel value of the histogram that corresponds to the clipping fraction becomes aligned 
with Xiow The judgement of clipping is done channel by channel, without regard to the 

25 amount of clipping in other color channels. Similarly, in the highlight region each of the 
control points Rhigh, Ghigh, Bhigh'that lies below Xhi g h is moved towards the reference point 
Xhigh, subject to the condition that the clipping fraction is not exceeded. Thus, if a 
particular channel does not exceed the predetermined clipping fraction its control point is 
aligned with Xhi g h. Otherwise, the channel value of the histogram that corresponds to the 

30 clipping fraction becomes aligned with Xhigh- Again the judgement of clipping is done 
channel by channel, without regard to the amount of clipping in other color channels. The 
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old and new positions of the control points are used to establish three linear look-up 
tables for correcting the color balance. 

Second Embodiment 

5 This embodiment, referred to as the color cast correction, is suitable for processing the 
majority of normal images that require color balancing and in addition is useful for 
images that exhibit strong color errors in the form of a distinct color cast. It can also be 
applied to the automatic correction of images without user intervention. The procedure 
will be described in two steps for clarity but, as will be appreciated by those skilled in the 

10 art, it can also be accomplished in a single step. In step 1 , for the shadow region each of 
the control points Ri ow , Gi ow , Bi ow is moved into exact alignment with the reference point 
Xiow (Note that one color control point is already coincident with Xi ow as it was used to 
determine Xi ow ). Similarly, in the highlight region each of the control points R h i g h, G h i g h, 
Bhigh is moved into exact alignment with the reference point X h i g h. This may result in one 

15 or more of the color channels being clipped by more than the permitted clipping fraction 
at one end or the other or both. In step 2, the three channel histograms are shrunk 
simultaneously by the same amount until none of the channels exhibits more clipping 
than the predetermined clipping fraction. The alignment of the three control points at 
each end of the histogram is maintained during this process. In other words, the Xi ow 

20 reference point is shifted upwards in value until no channel has more than the 

predetermined clipping fraction at the low end. Similarly, the X high reference point is 
shifted downward in value until no channel has more than the predetermined clipping at 
the high end. The old and new positions of the control points are used to establish three 
linear look-up tables, that is one for each of the colors, for correcting the color balance. 

25 

Third embodiment 

Following the correction according to the first or the second embodiment of the invention 
the color of the resulting image is well balanced. However, the true color of the scene, 
especially of neutral gray or white in the scene, may not correspond to the light source 
30 used to image the scene. Thus, a scene illuminated with summer skylight of color 
temperature 9500K would lead to a bluish white. The same scene illuminated with 
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average summer sunlight with some skylight at a color temperature of 6500 K would give 
rise to a neutral white. Illumination with a flash lamp of color temperature 5000K would 
produce a yellowish white and an incandescent lamp of 2900K would yield a more 
orange white. In some cases it is advantageous to restore the appearance of the original 
5 scene illuminant even though, strictly speaking, the colors are no longer balanced. There 
are many ways of achieving this objective as will be obvious to those skilled in the art. 
One niethod that has been found effective and convenient is the following. A control 
point is chosen for performing illuminant correction expressed, for example, as a 
lightness Lk in the CIE L*a*b* color space. While other opponent color spaces with a 

10 lightness or brightness axis may be used, color spaces such as CIE L*a*b* or CIE 
L*u*v* are preferred since their lightness representation is well matched to human 
vision. Though the exact value of this control point Lk is not critical provided it is not at 
the extremes of the lightness range, it is preferred to have a value near the middle of the 
lightness range for it is here that colors tend to have the highest saturation and 

15 consequently the effect of illuminant color correction is most easily discerned. An 

especially preferred value of Lk is that which corresponds to a mid-grey in the RGB color 
space, such as that represented by a neutral gray of R = G = B = 127, where the range of 
R, G and B lies between 0 and 255. A table is prepared of the CIE L*a*b* color a^, bk 
values for a black body radiator of this lightness (L k ) but of varying color temperature, k. 

20 Advantageously this table includes color temperatures from about 2000K to about 

10000K. It is sufficient to have about 5 to 15 values in the table, with more values at the 
lower temperatures, and to achieve other values by interpolation. Several interpolation 
methods may be used, but linear interpolation has been found both fast and adequate. The 
user selects an illuminant temperature and the corresponding a k , bk color components are 

25 retrieved. The shift in red, green and blue required to transform the CIE L*a*b* color 

with the coordinates (Lk, 0, 0) to that with the coordinates (Lk, a k , b k ) is then computed. A # 
gamma function correction is calculated accordingly for each channel using this adjusted 
control point and the lower and upper limits of the channel values. The correction 
function is expressed as a look-up table. The correction for illuminant is an optional one. 

30 Without this correction the image is correct for a light source of about 6500K color 
temperature. 
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Fourth embodiment 

Following correction according to one or other of the first three embodiments the image 
will have an appropriate color balance. However, the contrast of the image may have 

5 changed during color correction. While any contrast change may have led to a contrast 
improvement it is often desirable to modify only the color balance without changing 
contrast. This is especially desirable if the image contrast has already been independently 
adjusted. Accordingly, it is preferred to optionally correct the image contrast in such a 
way as to restore the original image contrast. This allows color to be balanced 

10 independently of other changes in the properties of the image. Many ways of 

accomplishing these objectives will occur to those skilled in the art. It has been found that 
the following procedure is both simple and gives satisfactory results. The lightness or 
brightness component of a suitable color space is computed for the original uncorrected 
image. Similarly the same lightness component is calculated for the image that results 

15 from color correction. A look-up table is constructed that maps the result image lightness 
histogram to the original lightness histogram, so restoring the original image contrast. 
While several color spaces, such as CIE L*a*b*, CIE L*u*v*, HSL, HVC or HSB, have 
a suitable lightness component for performing the correction, color spaces such as YIQ, 
YUV or YES are preferred because the transformation to and from these spaces is rapid 

20 and the Y brightness component is a good match to human perception. Of these, YIQ is 
especially preferred. 
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